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Abstract 

An explicit two-soliton solution for the derivative nonlinear Schrodinger equation 
with nonvanishing boundary conditions is derived, demonstrating details of inter- 
actions between two bright solitons, two dark solitons, as well as one bright soliton 
and one dark soliton. Shifts of soliton positions due to collisions are analytically ob- 
tained, which are irrespective of the bright or dark characters of the participating 
solitons. 
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The derivative nonlinear Schrodinger (DNLS) equation is an integrable model 
describing various nonlinear waves such as nonlinear Alfven waves in space 
plasma(see, e.g., [1,2,3,4,5,6,7]), sub-picosecond pulses in single mode optical 
fibers(see, e.g., [8,9,10,11,12]), and weak nonlinear electromagnetic waves in 
ferromagnetic [13], antiferromagnetic[14], and dielectric[15] systems under ex- 
ternal magnetic fields. Both of vanishing boundary conditions (VBC) and non- 
vanishing boundary conditions (NVBC) for the DNLS equation are physically 
significant. For problems of nonlinear Alfven waves, weak nonlinear electro- 
magnetic waves in magnetic and dielectric media, waves propagating strictly 
parallel to the ambient magnetic fields are modelled by VBC while those 
oblique waves are modelled by NVBC. In optical fibers, pulses under bright 
background waves are modelled by NVBC. 
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It has been known that integrable systems admit soliton solutions which pass 
through each other in a completely elastic fashion. The only traces of their 
interactions are shifts of positions and phases. In some systems, besides "pure" 
solitons which keep their forms unchanged in transmissions, there exist bound 
states of these pure solitons called breathers which periodically oscillate (see, 
e.g., [16]). 

Soliton solutions for the DNLS equation with VBC, including one-soliton 
solution[17] and multi-soliton formulas(e.g.,[18,19]), have been known. Re- 
searches on the DNLS equation with NVBC showed that its general one- 
soliton solution (corresponding to a complex discrete spectral parameter) is 
a breather which degenerates to a pure bright or dark soliton when the dis- 
crete spectral parameter becomes purely imaginary[20,21,4,5,22]. In known 
(l+l)-dimensional one-component integrable systems, the DNLS equation 
with NVBC is a rare instance which simultaneously supports pure bright soli- 
tons, pure dark solitons, as well as their bound states. Collisions between 
these solitons are thus important topics. However, like other NVBC prob- 
lems, a double-valued function of the spectral parameter inevitably appears 
in the inverse scattering transform (1ST) for the DNLS equation with NVBC, 
greatly complicating the 1ST. Early 1ST for the system performed on Riemann 
sheets only obtained modulus of a one-soliton solution [20] and asymptotic be- 
haviors of the modulus of an implicit pure two-soliton solution[21]. Although 
the phase of the one-soliton solution was found later, yielding a very com- 
plicated solution[4,5], it is a too tedious task to get an explicit multi-soliton 
solution based on the 1ST performed on Riemann sheets [20,21]. A recent 
multi-soliton formula using Backlund transformation [23] was also unable to 
explicitly demonstrate collisions between solitons. 

It has been suggested that constructing Riemann sheets for such NVBC prob- 
lems can be avoided if one performs the 1ST on the plane of an appropriate 
affine parameter [24]. The technique was recently applied to the DNLS equa- 
tion with NVBC, yielding not only a much simpler one-soliton solution than 
those in the literature but also a simple 1ST for further researches [22]. Immedi- 
ately following Ref . [22] , infinite number of conservation laws was derived by a 
simple standard procedure [25] and the evolution of a rectangular initial pulses 
in the system was considered, which was shown to be highly nontrivial and 
significantly different from all known results of other integrable systems[26]. 

It should be emphasized that these DNLS solitons might not exist in all of 
the physical systems mentioned above when the (l+l)-dimensional DNLS 
equation is not rigorously valid. For example, in space plasma, when the DNLS 
equation was generalized to a multi-dimensional one[27,28], it was shown that 
the dark solitons are not stable in presence of transverse perturbations [28]. 

In this letter, we only consider pure solitons of the DNLS equation with NVBC. 
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We modify the 1ST in Ref. [22] to the case when all discrete spectral parameters 
are purely imaginary and derive an explicit pure two-soliton solution. The 
solution consists of two bright solitons, two dark solitons, or one bright soliton 
and one dark soliton. Shifts of soliton positions due to collisions between them 
are obtained. We find that the collision between a bright soliton and a dark 
soliton is similar to usual collisions between two bright solitons or two dark 
solitons: the position of the faster one get a forward shift while the position 
of the slower one get a backward shift. We also show that the shifts of soliton 
positions are irrespective of the bright or dark characters of those participating 
solitons. 

We write the DNLS equation as 

iu t + u xx + i(\u\ 2 u) x = 0, (1) 

where the subscript denotes partial derivative. Its first Lax equation[17] is 
d x F = LF, (2) 




(3) 
(4) 

Here <7j(i = 1, 2, 3) are Pauli matrices, 

ffi=( 01 |, <x 2 = ° l> *s=| I (5) 



the bar stands for complex conjugate, and A is the time-independent spec- 
tral parameter. As there is no phase shift across the DNLS solitons with 
NVBC[4,22], the NVBC can be simply written as 

u — > p, as x — > ±oo, (6) 

where p is real. Asymptotic solutions of Eq. (2) is 

E ± (x, k) = (I- ipk- 1 a 1 )e- iXCxa3 , as x -> ±oo, (7) 
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where we introduce an affine parameter k to make ( = (A 2 + p 2 )" 2 , a double- 
valued function of A, become single- valued function of k, with 

\ = \(k-p 2 k- 1 ), C = \(k + P 2 k- 1 ). (8) 

We define Jost solutions, 

ty(x, k) — > E + (x, k), as x -> oo, (9) 
fc) — > fc), as x — > — oo, (10) 



where 

*(x,k) = (i/>(x,k), i/>(x,kj), 
$(x,k) = (<l>(x,k), 4>(x,k)), 

and the scattering coefficients by 

4>(x, k) = a(k)if)(x, k) + b(k)ip(x, k), 
<f>(x, k) = d(k)ip(x, k) — b(k)tp(x, k). 



(11) 
(12) 

(13) 
(14) 



tp(x,k), (f>(x,k) and a(k) are analytic in the first and the third quadrants of 
the complex k plane, while ip(x, k), <p(x, k) and a(k) are analytic in the second 
and the fourth quadrants. 

As shown in Ref. [22], on the plane of the affine parameter k, if k n \ = k n is a 
simple zero of a(k) in the first quadrant, then k n 2 = —k n , k n z = p 2 k~ l , and 
k n 4 = —p 2 k~ x are also simple zeros. For the case when all discrete parameters 
X n (n — 1, 2, . . . , N) are purely imaginary, all zeros of a(k) locate on the circle 
of radius p centered at the origin(p-circle), that is 

k n = pexp(i(3 n ), 0<f3 n <ir/2, n = 1, 2, . . . , N. (15) 

Then, k n3 = k nl , k n4 = k n2 , contributions of k n3 and k n ± must be dropped from 
relevant equations obtained in Ref. [22]. Therefore, for reflect ionless potentials, 
we get 

N 1.2 _ 1.2 

«(*) = ^ n (16) 

n=l K h n 
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where 77 = -2^ n/ 3 n , the inverse scattering equation, 



ip(x, k)e^ x = 



-trf 



N 



( 

I ' 1 

+ 2 E 

^ — ipk 1 e' iri+ ) n=i 



k n \ c n ip(x, k n ) 



k 



k 2 - k 2 



and the expression for reflectionless potentials (solitons), 

N 

u(x) = pe~ l2v+ - 2pe~ iv+ ^ fc n )e a " c " 3: . 

n=l ^ n 



Here 



X n = ip sin /3 n , ( n = p cos /3 n , 
u n = p 2 sin(2/3 n ), w n = p 2 (l + 2sin 2 /3 ra ), 



c n (i) = c n (0)e l2A " c "( 2A "-^ 2 )* = c n (0)e 



and 

1 °° 

Symmetric relations such as 

ip(x, p 2 k~ r ) = ip~ 1 kas^p(x, k), 

found in Ref. [22] are still valid. At k — k n , we have 
ip(x, k n ) = ip~ 1 k n a^(x, k n ), 

then, at k — k m , the first component of Eq.(17) is 



N 

ie ipm ^{x, k m ) = e -^ + e iXm ^ x + 2 ]T 



knCnlpljX-, k n ) j( An 



J 1.2 _ 1.2 
1 

It can be shown that for k n located on the p-circle, 
c n = ie ll3n x (a real number). 



We can set 

C n (0) = l Xn psHWnY^e V " X \ Xn = ±1. 



We also define 



n = is n (x-x n -v n t). (26) 



In principle, one can find ipi(x,k m ) by solving linear equations Eq.(23) and 
then get a multi-soliton solution with Eq.(20) and Eq.(18). We only consider 
one-soliton and two-soliton solutions in this letter. 

For the case of N = 1, we get 

Ul =pe-^t^ (27) 



where 

D 1 = 1 -i X ie il3l e- e \ (28) 

A 1 = 1 -i X ie- lWl e- e \ (29) 

1 °° D 
Vt = ^J(P 2 ~ \^i\ 2 )dx = tin (30) 

X 

and thus the one-soliton solution, 



ui = P-jjr = ( 31 ) 



which is identical to that obtained in the literature. It is a bright soliton 
for Xi = -1 or a dark soliton for Xi — 1- There is only one parameter, 
characterizing the soliton which is usually called one-parameter soliton [5]. 

For the case of iV = 2, we get 

u 2 =pe-^^, (32) 



where 



D 2 = 1 - i X ie i(3l e- ei - ix 2 e il32 e- d2 

-X ^K^l^U ^e- 9 ^, (33) 
sm (f3i + fa) 
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A 2 = l- ix\e~ iWl e' ei - i X 2^ iW2 e^ 2 

-XiX 2 Si 1^:^ e-^^)e^^. (34) 
sin ^ + f3 2 ) 



In order to find i] 2 , we find 

MD,)*^! - Im(Z) 2 )^MM = 6|D 2 p - |A 2 ft. (35) 
ax ax 4 



With this relation, we get 



% = 2 



f(p 2 - \u 2 \ 2 )dx = i In j^-, (36) 



a; 



and the two-soliton solution, 

U2=P^yT- (37) 



Assume (3 2 > Pi- At times long before collision (t — > — oo), in the vicinity of 
#i pa 0, #2 — > oo, w 2 ~ ui(^i), while in the vicinity of 2 ~ 0, 9i — > — oo, 
■u 2 ~ «i (6*2 + A), that is, 

« 2 «Ui(0i)+Ui(0 2 + A), (38) 



where 

A = 2 In 



sin(/3i + /3 2 ) 



sin(/3i - /? 2 ) 



> 0. 



(39) 



■u 2 consists of two well separated solitons, moving to the positive direction of 
the x-axis, with the slower soliton of Pi moving on the front. 

At times long after collision (t — > oo), in the vicinity of Q\ ~ 0, 9 2 — > — oo, 
■u 2 ~ -ui (0i + A), while in the vicinity of 6 2 ^ 0, Q\ — > oo, w 2 ~ «i(0 2 ), that is, 

m 2 ~ Mi(0i + A) + iii(0 2 ). (40) 



Now u 2 consists of two well separated solitons, with the faster soliton of p 2 
moving ahead. These asymptotic solutions show that the faster soliton of p 2 
catches up the soliton of P±, collides with the latter and leaves it behind, gets 
a forward shift Ax 2 = A/u 2 in position while the soliton of Pi gets a backward 
shift Axi = —A/ui in position. These shifts in position are independent of 
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Xi and X2, be., the bright or dark characters of the two participating soli- 
tons. They are in agreement with those in Ref.[21] obtained by discussing the 
asymptotic behavior of the modulus of an implicit two-soliton solution. 

Collisions between two bright solitons(Fig.l), two dark solitons (Fig. 2), and 
one bright soliton and one dark soliton(Fig.3 and Fig.4) are graphically shown, 
where shifts of soliton positions due to collisions are vividly seen. We choose 
X\ = A/(2vi) and x 2 = A/(2z/ 2 ) in these figures so that the two solitons 
completely overlap at (x,t) = (0,0). 

In summary, we find an explicit two-soliton solution for the DNLS equation 
with NVBC. Shifts of soliton positions due to collisions between solitons are 
analytically obtained. Details of collisions between two bright solitons, two 
dark solitons and one bright soliton and one dark soliton are graphically shown. 
It is interesting to note that these shifts only depend on parameters of the 
participating solitons, irrespective of their bright or dark characters. Observing 
Eq.(30) and Eq.(36), one can find that the relations between rj and D for one- 
soliton and two-soliton are the same. The relation can possibly be extended 
to multi-soliton case and be helpful to find an explicit multi-soliton solution. 
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Fig. 1. Collision between two bright solitons, p = 1, f3i = 7r/12, = 7r/24. Variables 
in the figure are dimensionless. 




Fig. 2. Collision between two dark solitons, p = 1, j3\ = 7r/15, = tt/6. Variables 
in the figure are dimensionless. 
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Fig. 3. A bright soliton of fa = 2it/15 chases a dark soliton of fa = 77r/60, p = 1. 
Variables in the figure are dimensionless. 




Fig. 4. A dark soliton of fa = tt/6 chases a bright soliton of fa = tv/12, p = 1. 
Variables in the figure are dimensionless. 
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